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We demonstrate and characterize a high-flux beam source for cold, slow atoms or molecules. The 
desired species is vaporized using laser ablation, then cooled by thermalization in a cryogenic cell 
of buffer gas. The beam is formed by particles exiting a hole in the buffer gas cell. We characterize 
the properties of the beam (flux, forward velocity, temperature) for both an atom (Na) and a 
molecule (PbO) under varying buffer gas density, and discuss conditions for optimizing these beam 
parameters. Our source compares favorably to existing techniques of beam formation, for a variety 
of applications. 

PACS numbers: 39.10.+J 
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Cold, slow beams of atoms or molecules are of wide 
utility. A common use of such beams is as a source 
for loading into traps, where the particles can be fur- 
ther cooled and manipulated, e.g. to create Bose con- 
densates or Fermi degenerate gases. Because the number 
of trapped particles is typically limited by the character- 
istics of the initial beam (flux, forward velocity, tempera- 
ture, etc.), significant effort has been put into developing 
improved atomic beam sources 0, 0, IE & IE IE ■ Work 
on developing cold molecular sources has recently been a 
particularly active field of research [1 IE EE El El El 
As with atoms, one of the aims is to produce quantum 
degenerate gases, including those comprising strongly in- 
teracting electric dipoles [fl 0, 0] . 

For the purposes of loading traps, an ideal source would 
produce a large flux of any atom or molecule at temper- 
atures less than the depth of the trap, Tj. For most 
currently used trap technologies, T f < 1 K. For species 
amenable to laser cooling, this temperature is within the 
capture range of a MOT E3- F° r paramagnetic atoms 
and molecules, strong magnetic traps can have depths ex- 
ceeding IK |l8|. For polar molecules, electric field-based 
traps can reach similar depths |l9|, |2CJ . 

We describe here a simple, robust source that can op- 
erate with nearly any atomic or molecular species, and 
which produces a beam at high flux with translational 
and rotational temperatures near 1 K. We demonstrate 
this source for both an atom (Na) and a polar molecule 
(PbO), and discuss the different regimes of beam for- 
mation. Our analysis makes it possible to estimate and 
optimize various characteristics of the beam source for 
general use. We believe this provides an attractive alter- 
native to many beam techniques now in use. 

A simple outline of the operation of our source is 
as follows. Atoms or molecules of the desired species, 



A, are first vaporized by laser ablation. This produces 
N particles of A per pulse, at a high temperature Ti. 
The ablation takes place inside a cryogenic cell filled 
with He buffer gas at low temperature Xt> and density 
riHe- After a characteristic number of collisions Af, 
the translational temperature T of A comes arbitrar- 
ily close to equilibrium with the buffer gas, such that 
T « (1 + e)T b when N = ~Kln{eT b /T t ) 0. [Here 
k = [vtiA + rriHe) 2 1 (ZrnAmHe), and m# e (m^) is the 
mass of He (A)]. Rotational degrees of freedom are also 
cooled during these collisions. By allowing both He and 
A particles to exit the cell via a small hole, a beam of 
A is formed. The beam persists for a duration given by 
the diffusion lifetime of A in the cell, which is limited by 
sticking of A particles upon contact with the cell walls. 

The number of cold particles of A in the beam is de- 
termined by both riHe and the cell geometry. During 
thermalization, a particle of A typically travels a dis- 
tance R ~ M I (jiHe^t), where a t is a thermally-averaged 
cross section for elastic collisions. Hence, for a cell with 
distance Rh from the ablation point to the hole, the par- 
ticles of A will be efficiently thermalized before exiting 
the hole only if R < Rh . In addition, the purely geomet- 
ric probability for a particle of A to escape in the beam 
is governed by the ratio d/ Rh, where d is the diameter of 
the exit hole. 

The forward velocity, Vf, of the thermalized beam of 
A particles also depends on both n# e and d. Specifically, 
Vf is determined by the ratio d/X, where A = l/(n# e cr c ) 
is the mean free path of A particles in the cell; here a c is 
the elastic cross-section for cold A-He collisions. In the 
effusive limit (A > d), u/ will be given approximately by 
the thermal velocity of cold A particles, i.e. Vf = 
y / 2kBTb/mA- By contrast, when A « the A particles 
will become entrained in the outward flow of He, so that 
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FIG. 1: A schematic of the beam apparatus. When detecting 
fluorescence, a lens (not shown) collimates a fraction of the 
fluorescence light and directs it out a window to a PMT. 



Vf rs VHe = y/ZkBTb/mjie- Since vha S> mue for most 
species of interest, Vf is much smaller in the effusive limit 
than for an entrained beam. 

Note that the conditions for efficient thermalization 
and for a slow beam are in conflict. Thermalization is 
most efficient for riHe above a threshold value, but ef- 
fusive flow demands that nue be less than a typically 
different threshold. The highest flux of cold, slow A par- 
ticles is obtained when Rh and d are chosen so that these 
thresholds in buffer-gas density coincide (namely, when 
Rh/ d k, Afa c /at), and n# e is set at this common thresh- 
old value (njj], = da c = Rh&t)- 

A schematic of the beam setup is given in Fig. ^ 
The buffer gas cell is a brass box ~10 cm on edge. The 
cell is mounted in vacuum, with the top face attached 
to the cold plate of a liquid He cryostat. An 0.8 mm 
long exit hole with d = 3 mm is centered on one side 
face. The bottom and other side faces are covered with 
windows for optical access. Several ablation targets are 
mounted on the top face at Rh ~ 6 cm. The ablation 
light consists of laser pulses of ~ 5 ns duration, with 
energy ~ 15 mJ, focused to a spot size < 1 mm, at a 
wavelength of 532 nm. The ablation laser is typically 
fired at 10 Hz repetition rate. We produce Na atoms with 
sodium metal or NaCl targets, and PbO molecules with 
a vacuum hot-pressed PbO target. With our ablation 
conditions, typically Tj « 1000 K for both species. 

Buffer gas continuously flows into the cell through a 
narrow tube which is thermally anchored to the cold 
plate. This ensures that the buffer gas is at the tem- 
perature of the cold plate upon entry to the cell. The 
cell walls, and thus the He gas, are typically at T& = 5 K. 
Helium in the cell must be replenished as it flows con- 
tinuously through the exit hole. A calibrated metering 
valve at room temperature is used to control the flow 
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into the cell, and hence the density, of the buffer gas. 
We determine nue to within a factor of 2. Good vac- 
uum is maintained in the beam region by means of a 
coconut charcoal sorption pump with a pumping speed 
of -1000 Is -1 . 

For a typical elastic cross-section a c « 3 x 10 15 cm -3 , 
the crossover between effusive and entrained flow of 
A — i.e., the condition d = A — occurs for riHe — n c « 
10 15 cm -3 . This should be compared to the density 
= n t required for full thermalization of A parti- 
-i.e., such that R « N / '{nHe<^t) = Rh- Assum- 
ing cr t w cr c , we find that for Na in our cell, nt ~ n c . 
Thus it should be anticipated that our cell is near the 
optimal geometry for producing a maximal flux of slow, 
cold Na. By contrast, the larger mass of PbO makes 
J\f much larger than for Na, implying that our cell ge- 
ometry is not optimal for PbO. We characterize the the 
beam source for both species within a range of densities 
around the anticipated optimal condition for Na, namely 
n He w 0.2 - 5 x 10 15 cm" 3 . 

The beam source is monitored using laser spectroscopy. 
Doppler shifts and widths of the spectra are used to de- 
termine beam velocity profiles. Signal size and timing 
yield the column density and particle dynamics. To mea- 
sure the longitudinal (transverse) velocity profile, a probe 
laser beam is sent collinear with (perpendicular to) the 
molecular beam. For Na, a third probe beam monitors 
atoms inside the cell. 

We monitor Na atoms via absorption of a probe laser 
tuned to the 3si/ 2 —> 3p 3 / 2 transition at wavelength 
Ajva = 589 nm. The probe laser frequency is continu- 
ously scanned over the entire absorption profile (a range 
of ~ 1 GHz), at a rate of 1 kHz, resulting in several com- 
plete spectral profiles of the Na atomic beam for each 
ablation pulse. The relatively slow scan of the laser en- 
sures that we are sensitive only to nearly-thermalized Na 
atoms; hot atoms traverse the detection region before a 
scan is complete. 

PbO is monitored via laser-induced fluorescence. This 
probe laser is tuned to the X(v" = 0) — >B(t/ = 5) transi- 
tion at \pbo,e — 406 nm. Fluorescence is detected using 
a photomultiplier tube (PMT), with interference and col- 
ored glass filters to selectively observe the B(t/ = 5) — * 
X(v" = 4) transition at Xpbo.f = 460 nm. A shutter with 
an opening time of — 7 ms is used to shield the PMT from 
the initial glow following each ablation pulse. Again, 
this technique ensures that only slow-moving molecules 
are detected. The signal is averaged over several (typi- 
cally 10) shots of the ablation laser with the probe laser 
frequency fixed. A single spectral scan consists of fre- 
quency points separated by 10-30 MHz across a span of 
up to several GHz. Our quantitative determination of 
molecule number from fluorescence counts includes esti- 
mates of the detection efficiency and branching fraction 
for the detected transition |22|. We estimate the error 
in our absolute number measurement of PbO to be less 
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FIG. 2: The number of cold A particles, Nh, that emerge in 
the beam as a function of riHe ■ Curves with specific functional 
forms have been inserted to show the different scaling regimes. 
The intersection of the curves corresponds to the condition 
R — Rh, where a typical A particle is thermalized just as it 
reaches the hole. 

than a factor of 2. 

For Na, the in-cell probe beam was used to determine 
both N and a c . We find N Na ks 10 14 /pulse for both 
the Na and NaCl targets. We measure diffusion lifetimes 
of r [ms] ~ 4 x 10 -15 x riHe [cm -3 ]. From this we in- 
fer <7 C rs 3 x 1CP 15 cm 2 . For PbO, previous work has 
measured an ablation yield of w 10 12 /pulse [2^, and our 
measurements with an in-cell probe indicate a compara- 
ble yield. 

In Fig. |2we plot the number, Nh, of thermalized par- 
ticles of species A exiting the hole as a function of nn e - 
We also show the results of a Monte Carlo simulation of 
the beam formation process. The condition for full ther- 
malization is apparent in both the experimental and sim- 
ulated data for Na, while as expected we do not appear 
to reach the condition of full thermalization for PbO. 

In the Na data and simulations, we find that Nh in- 
creases rapidly (approximately oc n He ) up to a critical 
value of riHe, above which Nh is roughly constant at 
its maximum value Nh, max - The low-density scaling is 
consistent with a simple picture in which A particles 
are distributed uniformly over a volume of characteris- 
tic length L t ~ MX, when they have thermalized to near 
Tb- (This broad distribution arises, in our simulations 
of the thermalization process, from the spreads both in 
the number of collisions required to thermalize and in 
the free path between collisions.) At high density, the 
condition Nh — Nh.max arises in our simulations when 
R < Rh for essentially all A particles. In this regime, 
the fraction of A particles escaping, f max = Nh, m ax/N, 
is given roughly by the area of the hole to the area of a 
hemisphere at radius Rh, i.e., f max ~ <i 2 /(8i? 2 ). In our 
geometry, f max rj 3 x 1(T 4 . 

The simulated beam data in Fig. 121 is scaled to match 
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FIG. 3: PbO and Na beam mean forward velocities, Vf, as a 
function of nHe- Extrapolation of the data to zero buffer gas 
density is illustrated by best-fit lines (dashed). 

the experimental Na data by adjusting the values of iVjva 
and <T t - The resulting value, N^ a rs 5 x 10 13 , is in 
reasonable agreement with the determination from the 
in-cell probe. The fitted thermalization cross-section, 
(T(~lx 10 -15 cm 2 , is somewhat smaller than a c . This 
is reasonable, since elastic cross sections are typically 
smaller at higher collision energies [24^. 

Fig. |31 shows the average forward velocity, v/, of the 
beams of A particles as nue is varied. For both Na and 
PbO, the data show a nearly linear increase of Vf with 
riHe, with the velocity of the lighter species always larger. 
This behavior is consistent with the following simple pic- 
ture. A slowly-moving particle of A takes a time T e to 
exit the hole, where T e ~ d/vA- During this time, it un- 
dergoes N e collisions with fast, primarily forward- moving 
He atoms, where N e ~ n^ e cr c w^ e T e . Each collision im- 
parts a momentum transfer ApA ~ mHe v He- This results 
in a net velocity boost Ava, given by Ava ~ VAd/X oc 
nHe- This picture should be roughly valid for densities 
below the regime of full entrainment, where Vf ~ vn e - 
The velocities we measure for Na are approximately re- 
produced by modeling of the beam formation process 
with our measured value of a c . 

This picture also predicts that the behavior of Vf, when 
extrapolated to njje = 0, should yield the velocity of an 
effusive beam of A particles at temperature Tf,. To make 
this comparison, it is critical to note that our detection 
technique is sensitive to molecules within a roughly cylin- 
drical volume, of diameter Dd and with length L extend- 
ing from the exit hole. Under our conditions, where Dd ~ 
d < L, it can be shown that an effusive beam will exhibit 
a velocity distribution close to f(v) oc v 2 e ~ m AV /{2k B T h ) 
and a mean velocity v e ff ~ I.ISva- Our extrapolated 
data is within ~25% of this prediction for both species. 

Fig. 0] shows the temperature of the beams vs. % e , 
as measured by the velocity spreads in the longitudi- 
nal and transverse directions, as well as the rotational 
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FIG. 4: PbO and Na temperatures as a function of tlhh as 
determined by fits to spectral profiles. "Transverse", "longi- 
tudinal" , and "in cell" refer to the location of the probe laser 
(see Fig. 0. Typical errors for the transverse and in cell Na 
temperatures are comparable to the marker size. 



population distribution (for PbO). For transverse tem- 
perature we fit to a distribution of the form ft(v) oc 
e -m A v /(2fc B T) -p Qr i on gitudinal temperature we use 
fl(v) oc e -^A(v^v l f/(2k B T) ^ Fits of the i on gi tuc ii na i d ata 
to the effusive distribution fi, e ff{v) oc v 2 e ~ m = v /(2*bT) 
were poor, consistent with the partial cntrainmcnt of A 
in the helium flow. For the rotational temperature, we 
use Clebsch-Gordan and Honl-London factors to deter- 
mine the relationship between fluorescence intensity and 
initial state population for various rotational lines. The 
ratio of initial state populations for 2-3 rotational lines 
determine the temperature corresponding to a Boltz- 
mann distribution. Our data indicates complete ther- 
malization of all detected particles. Note that we observe 
no additional cooling below Tt>, as would be expected for 
entrainment in a fully supersonic He flow. 

The beam source described here can be readily adapted 
to the needs of a wide range of experiments. For exam- 
ple, it can be used for a wide variety of species with 
performance similar to that described here. The total 
flux depends linearly on the ablation yield N, which for 
any given target is difficult to predict a priori. How- 
ever, for nearly every species we have tried (here and 
in many related experiments), it has been possible to 
achieve large values of N by a suitable choice of precursor 
material. For example, under operating conditions simi- 
lar to those used here, we have obtained N = 10 12 — 10 14 
for a variety of metal atoms and iV = 10 11 — 10 13 for 
many s pec ies of diatomic molecules (including radicals) 
|25U2q . l27ll2S| . In addition-subject to limitations of cool- 
ing power or gas load-the source could be run at lower 
temperatures, higher repetition rates, or with higher ex- 
traction efficiency (e.g., by using several separated exit 
holes) . Our discussion of the beam formation mechanism 



makes it straightforward to determine the effects of such 
changes. 

A key possible improvement to the source would be 
the addition of a guide, either magnetic (for paramag- 
netic species) or electric (for polar molecules). In both 
cases, He is unaffected by the guide potential and will 
exit through the sides of the guide, allowing extraction 
of the A beam into a region of ultra-high vacuum. Under 
the conditions described here, this beam source could be 
used to load a peak flux of up to 2 x I0 11 s _1 Na atoms 
into a simple permanent magnet guide such as the type 
described in Ref. hf. An electrostatic guide such as that 
described in Ref. |l2( could be loaded with a peak flux of 
~ 10 9 s _1 PbO molecules in the J=l rotational state, ad- 
equate for loading into a microwave |l9) or electrostatic 
[2jj trap. The high fluxes from our source could result in 
substantial improvements in atomic and molecular trap- 
ping experiments that depend on large initial numbers. 

We acknowledge support from the National Science 
Foundation under Grant Nos. DMR-0325580 and PHY- 
0071311, the Army Research Office, the W.M. Keck 
Foundation, and the David and Lucile Packard Founda- 
tion. 



[1 

[2. 

[3 
[4 

[5. 
[6 
[7; 
[8 

[9 

[10 

[11 

[12 
[13 

[14 
[IS! 

[16 
[17 
[18 

[19 

[20 
[21 

[22] 



and 



105, 



Z. T. Lu et al., Phys. Rev. Lett. 77, 3331 (1996). 

C. Slowe, L. Vernac, and L. V. Hau (2004), 
physics/0407040. 

J. Schoser et al., Phys. Rev. A 66, 023410 (2002). 
K. Dieckmann, R. J. C. Spreeuw, M. Weidemuller, 
J. T. M. Walraven, Phys. Rev. A 58, 3891 (1998). 
P. Cren et al., Eur. Phys. J. D 20, 107 (2002). 
R. S. Conroy et al, Opt. Commun. 226 (2003). 
E. Nikitin et al., Phys. Rev. A 68, 023403 (2003). 
M. Gupta and D. Herschbach, J. Phys. Chem. A 
1626 (2001). 

H. L. Bethlem, G. Berden, and G. Meijer, 
Phys. Rev. Lett. 83, 1558 (1999). 

M. R. Tarbutt et al., Phys. Rev. Lett. 92, 173002 (2004). 
M. S. Elioff, J. J. Valentini, and D. W. Chandler, Science 
302, 1940 (2003). 

S. A. Rangwala et al., Phys. Rev. A 67, 043406 (2003). 
R. Fulton, A. I. Bishop, and P. F. Barker, 
Phys. Rev. Lett. 93, 243004 (2004). 

L. Santos, G. V. Shlyapnikov, P. Zoller, and M. Lewen- 

stein, Phys. Rev. Lett. 85, 1791 (2000). 

M. A. Baranov, M. S. Mar'enko, Val. S. Rychkov, and 

G. V. Shlyapnikov, Phys. Rev. A 66, 013606 (2002). 
M. Baranov et al., Phys. Scr. 102, 74 (2002). 

E. L. Raab et al., Phys. Rev. Lett. 59, 2631 (1987). 

J. G. E. Harris et al., Review of Scientific Instruments 

75, 17 (2004). 

D. DeMille, D. Glenn, and J. Petricka, Eur. Phys. J. B 
31, 375 (2004). 

H. Bethlem et al., Nature 406, 491 (2000). 

R. deCarvalho et al., European Physical Journal D 7, 289 
(1999). 

E. A. Dorko et al., Chem. Phys. 102, 349 (1986). 



■5 



[23] D. Egorov et al., Phys. Rev. A 63, 030501(R) (2001). 

[24] L. Landau and E. Lifshitz, Quantum Mechanics (Non- 
relativistic Theory), vol. 3 of Course of Theoretical 
Physics (Pergamon Press, 1977), third revised ed. 

[25] J. Kim et al., Phys. Rev. Lett. 78, 3665 (1997). 

[26] B. Friedrich, J. Weinstein, R. deCarvalho, and J. Doyle, 



J. Chem. Phys. 110, 2376 (1999). 
[27] B. S. Zhao et al., Rev. Sci. Instrum. 75, 146 (2004). 
[28] C. I. Hancox et al., Nature 431, 281 (2004). 
[29] H. L. Bethlem et al., Phys. Rev. A 65, 053416 (2002). 



